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The appearance of travelling-wave-type solutions in pipe Poiseuille flow that are dis- 
connected from the basic parabolic profile is numerically studied in detail. Wc focus on 
solutions in the 2-fold azimuthally-periodic subspace because of their special stability 
properties, but relate our findings to other solutions as well. Using time-stepping, an 
adapted Krylov-Newton method and Arnoldi iteration for the computation and stabil- 
ity analysis of relative equilibria, and a robust pseudo-arclength continuation scheme we 
unfold a double-zero (Takens-Bogdanov) bifurcating scenario as a function of Reynolds 
number (Re) and wavenumber (n). This scenario is extended, by the inclusion of higher 
order terms in the normal form, to account for the appearance of supercritical modulated 
waves emanating from the upper branch of solutions at a degenerate Hopf bifurcation. 
These waves are expected to disappear in saddle-loop bifurcations upon collision with 
lower-branch solutions, thereby leaving stable upper-branch solutions whose subsequent 
secondary bifurcations could contribute to the formation of the phase space structures 
that are required for turbulent dynamics at higher Re. 



1. Introduction 

Subcritical transition and the sustainmcnt of turbulence at moderate flow rates in shear 



flows is a problem of great theoretical complexity and practical relevance (Grossmann 
2000| |Eckhardt et aL||2007| |Eckhardt||2009[ ). Since Osbourne [Reynolds| ( p83] ) much cele 



brated work on the onset of turbulent motion in a straight pipe of circular cross-section, 
this problem has become a benchmark of transition without linear instability, i.e. sub- 



critical transition. Many theoretical 


(Boberg & Brosa||1988 


Brosa & Grossmann||1999), 


numerical (Schmid & Henningson 


1994 


Zikanov 1996 Shan et al. 1999) and experi- 


mental (Wygnanski & Champagne 


1973 


Wygnanski et al. 


1975 Darbyshire & Mullin 


|1995| Hof et a/.|2006) efforts have been devoted at comprehending how the laminar state 



(Hagen-Poiseuille flow) becomes unstable to finite amplitude perturbations despite its 
linear stability ( Pfennigerj|1961 Meseguer fc Trefethen|[2003 1 . 



The numerical computation of finite amplitude secondary solutions in the form of trav- 



elling waves (Faisst & Eckhardt 2003 Wedin & Kerswell 2004), whose existence relies 



on a mechanism called the self-sustained process (Waleffe 1997), initially advanced to 
account for turbulent regeneration, has brought renewed interest to the problem within 
the dynamical systems community. These waves typically appear in saddle-node bifurca- 



tions and break some of the symmetries of the problem while retaining others (Pringle 



& Kerswell 2007 Pringle, Duguet & Kerswell 20091. Their lower branches have been 



repeatedly shown to dwell on the critical threshold (Duguet, Willis & Kerswell 2008 



Mellibovsky & Meseguer 2009 1 separating the basins of attraction of laminar and tur- 



bulent flows, which suggests their implication in the transition process. Upper-branch 
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solutions, with wall frictions closer to turbulence, seem to play a role in developed tur- 
bulence ( Schneider, Eckhardt fe Vollmer||20"07 Kerswell fc Tutty]|2007| ) and have been 
observed in experiments ( Hof et aZ.||2004 ). 

The 2-fold discrete azimuthally-periodic subspace constitutes a special case in that it 
contains a family of travelling waves with unusually low-dimensional unstable manifolds. 
All other known families of travelling waves exhibit at least one unstable mode within 
their subspace that is shared by all waves, thereby rendering the family globally unsta- 
ble. This is not the case of a specific shift-reflect, 2-fold azimuthally-periodic family of 
travelling waves, the special stability properties of which can be qualified as convenient. 
Lower-branch solutions of this family possess a single unstable direction within the sub- 
space, which ultimately assists bare symmetry-restricted time evolution in converging 
them numerically (Duguet, Willis & Kerswell 20081, when adequately combined with 
edge tracking techniques ( |Skufca, Yorke fe Eckhardt|2006 Schneider, Eckhardt & Yorke 



[2007| . The lower branch being a saddle, it would not be surprising if the upper branch 
constituted a stable node, at least for some of the parameter values. Additional unstable 
directions arise when analysing the stability in full space or when considering subhar- 
monic instability to longer wavelengths, but the convenient stability properties within 
their azimuthal subspace, which can by itself sustain turbulence, renders it a particularly 
nice playground for the search of all sorts of time-dependent solutions with simple nu- 
merical machinery. These solutions may help understand turbulence within the subspace 
and give hints as to what might be happening in full space. 

This study aims at analysing and describing the phenomena associated to the emer- 
gence of disconnected solutions in shear flows. While the current analysis focuses on a 
given family of travelling waves in circular pipe flow, the results are thought to be com- 
mon feature in a vast range of other travelling wave families in pipe flow and several 
equilibria and relative equilibria in a variety of other shear flows such as plane Couette 
(lNagata|1997[|Wang , Gibson fc Waleffe|2007[) or plan e Poiseuille ( |Pugh fc Saffman|1988 



Soibelman & Meiron||1991l lEhrenstein & Koch||1991|) flows. 



The outline of the paper is as follows. In f|2j we present the pipe Poiseuille flow prob- 
lem and sketch a numerical scheme for the integration of the resulting equations. The 
symmetries of the problem are discussed and numerical methods for the computation 
and stability analysis of relative equilibria are introduced. Some aspects of normal form 
reduction in the case of bifurcating relative equilibria are also reviewed in |2j The main 
results of the computations are presented in where the different types of solutions 
resulting from an extensive parameter space exploration are exhibited, and evidence of a 
double-zero bifurcation scenario is exposed. These results are thoroughly discussed and 
explained in Q by considering the inclusion of higher order terms in the relevant normal 
form. Conclusions and ensuing prospects and challenges are summarised in Sp] 



2. Formulation and methods 

2.1. Equations and numerical scheme 

The constant flow-rate motion of an incompressible Newtonian fluid through a straight 
pipe of circular cross-section is considered. The governing parameter is the Reynolds 
number, defined as Re = UD/v, where U is the mean axial flow speed, D the pipe 
diameter and v the kinematic viscosity of the fluid. The dynamics is governed by the 
incompressible Navier-Stokes equations. 

In cylindrical nondimensional coordinates x = (r,8,z), the basic Hagen-Poiseuille so- 
lution of the Navier-Stokes problem reads u b = r + 9 + z = (1 — r 2 ) z. The 
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Navier-Stokes equations for the velocity-pressure u = (u ri ug,u z )-p perturbation fields 
are 

' ' ' ! • V)(tl b + U) - (Mb • V)« + / Z, 



d t u 



-Vp 



Re 



-Am 



V • m = 0, 



(2.1) 
(2.2) 

(2.3) 
(2.4) 

where an adjustable axial forcing / = f(t) in (2.1) ensures the constant mass-flux con- 
straint (2.3). The non-slip boundary condition at the wall and the axial and azimuthal 



p2tt p 1 

Q(u) = (u-z) rdrdO = 0, 

JQ JO 



it(l, 8, z; t) = 0, u(r, 6 + 27r/n s , z, t) — u(r, 6,z + 2tt/k, t) — u(r, 8, z, t), 



periodicities are enforced by (2.4), where k is the non-dimensional wavenumber, here con- 



sidered as an additional parameter selecting the axial periodicity of the flow structures. 
In this study, the axial wavenumber is varied in the range k £ [1.1, 1.85], corresponding to 
pipe periodicities A £ [3.40, 5.71], in units of D/2. This range covers all relevant existing 
solutions up to Re ~ 1500 while discarding subharmonic instabilities that might compli- 
cate the bifurcation scenario under examination at higher Re. The azimuthal wavenum- 
ber is set to n s — 2 throughout, meaning that only 2-fold azimuthally-periodic fields 
are considered. All solutions living in this subspace are also solutions to the full Navier- 
Stokes equations with n s = 1, but the instabilities associated to perturbations without 
the n s = 2 symmetry are suppressed, rendering some of the solutions stable and thus 
accessible through time evolution. 

For the spatial discretisation of ( 2.1||2.4 ) we use a solenoidal spectral Petrov-Galerkin 
scheme thoroughly described and tested by Meseguer fc Mellibovsky| ( |2007[ ). The method 
implicitly fulfils divergence free and boundary conditions, and dispenses with the need 
for dealing with the pressure field. The velocity field is expanded following 



JV 



u(r,6,z;t)= E 



-i(K,lz-\-n s n&) 



ui n (r;t), 



(2.5) 



l = -Ln=-N 



M 



(1,2) 
nm 



r), 



(2.6) 



(1 2) 

with O; ' ' the complex expansion coefficients defining the vector state a. In this way, 
the field decouples into its axial-azimuthal Fourier components ui n . As already noted, 
the expansion basis vi nm (r), presented in the work of jMeseguer fc Mellibovsky] (|2007| eq. 



(24-27)), is chosen so that boundary conditions (2.4) and solenoidality (2.2) are implicitly 



satisfied. The spectral resolution, checked as adequate for the computations performed 
in this study, has been set to L = 16, N = 12 and M — 36, corresponding to ±16 
axial and ±24 azimuthal Fourier modes, and to 37 Chebyshev collocation points for the 
radial coordinate. Upon projection on the test basis Vi nm (r) (Meseguer & Mellibovsky 



2007 eq. (32-35)), we are left with a system of ordinary differential equations for the 
complex- valued state vector a: 



A a 



1 

lie' 



Si + 6 2 a b{a, a) + fC, 



(2.7) 
(2.8) 

where the subscripts in Meseguer & Mellibovsky (2007 eq. (47)) have been omitted for 



Q • a = 0, 



simplicity and B has been decomposed into two separate contributions, Bi corresponding 
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to the viscous dissipation term, with the dependence on Re rendered explicit, and 62 to 
the linear advection terms. C is the projection of z onto the test basis and is therefore 
responsible for keeping the mass-flux constant, represented by the linear scalar equation 



(2.8), by instantaneous adjustment of /. 



The system (2.7) is evolved in time using a 4th order linearly implicit method with 
At = 1 x 10~ 2 , where time is measured in units of D/(4U). 

It will be handy to define some global properties of the flow that are independent of 
solid-body motions such as translation or rotation. One such quantity is the normalised 
energy of any given perturbation field it, which is given by 



1 



2e b Jo 



2-jt/k, p2-K pi 



JO 



e(u) = — / dz d6 / rdr u* ■ u, (2.9) 



with eb = 7t 2 /(3k) the energy of the basic flow and * symbolising complex conjugation. 
This energy decouples exactly into the sum of its axial-azimuthal Fourier components: 

L N 

e(u)= ^ ^2 £ ln = £ 1D + £2t + £2z + £3D, (2.10) 
l=-Ln=-N 

where ei n = e(ui n ) is the energy associated to the (I, n) axial-azimuthal Fourier mode. To 



the right of (2.10) the Fourier modes have been grouped in meaningful sets corresponding 
to the axisymmetric streamwise- independent component (eiDj mean flow component), 
the non-axisymmetric streamwise-independent component (e2z, energy signature of vor- 
tices and streaks), the axisymmetric streamwise-dependent component (£20 > mean wavy 
component) and the non-axisymmetric streamwise-dependent component (£30, purely 
three-dimensional energy) . 

Another global property is the mean axial pressure gradient that is required to drive 
the flow at constant mass-flux, relative to its laminar flow value. It can be easily stated 
in terms of the adjusting intensity of the axial forcing / and of Re: 

(y p ), = i + M. ( 2.ii) 

Finally, it will be useful to define the instantaneous axial and azimuthal phase speeds 
for any given velocity field. Their definition is clear for solid-body moving solutions, but 
can be extended to any velocity field by decomposing the time-evolving state vector into 
two contributions: one associated to solid-body motion and the other to a modulation. 
This decomposition can be made unique by minimising, as a function of the phase speeds, 
at any given time, the 2-norm of the modulational component: 



mm 

Cz,Cg 



a(t + At) - o(t) e -^ Klc ' +n " nc ^ At (2.12) 



For all rotating and travelling waves, c z and eg are constant, the minimum (2.12) is 
exactly zero, and they coincide with the axial and azimuthal phase speeds, respec- 
tively. The effective computation is done instantaneously by considering two consecutive 
time instants separated by At and by applying a root-finding Newton method to the 



2-dimensional Jacobian of the norm in (2.12) 



2.2. Symmetries of the problem 

Pipe flow and its basic solution possess 0(2) x 50(2) symmetry, i.e. it is invariant under 
all azimuthal rotations about the axis and reflections with respect to all diametral planes, 
as well as under all axial translations. All travelling wave families discovered so far are pe- 
riodic and have axial and azimuthal dependency, thus breaking all continuous symmetries. 
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While, in the azimuthal direction, the natural periodicity of the pipe renders all solutions 
invariant under rotations of 2tt, travelling waves can have an additional discrete symme- 
try that renders them invariant under the cyclic group Z ns = {.fiW, #27i72i ^2?r/n s } 
with n s the azimuthal wave number and 

R a (u)(x) = R a (u, v, w)(r, 6*, z; t) — (it, v, w)(r, 9 + a, z;t). (2.13) 

In the axial direction, on top of the trivial shift by a wavelength, the axial symmetry rup- 
ture gives rise to a special type of time periodicity, in the sense that the time dependence 
is a bare drift in the direction of the SO{2) symmetry. Travelling waves are therefore 
best described as relative equilibria that posses the continuous space-time symmetry 

(u, v, w) tw (r, 9, z; t) = (u, v, w) tw (r, 9,z- c z t; 0), (2.14) 

where c z is the axial drift speed. In a comoving reference frame travelling downstream 
with speed c 2 , travelling waves appear as stationary solutions. Near a relative equilibrium 
the drift dynamics is trivial and decouples from the dynamics orthogonal to the drift. As 
a result, bifurcations of relative equilibria can be analysed in two steps, first describing 
the bifurcations associated to the orthogonal dynamics, then adding the corresponding 
drift along the travelling direction ( Krupa|fl990 1. 



There are a couple of additional non-trivial symmetries that leave some of the known 
travelling waves unaltered. The first one is a mirror symmetry with respect to n s diame- 
tral planes: 

Mu(x) = M(u,v,w)(r,6i + 6,z;t) = (u, -v, w)(r, 4 - 6, z; t), (2.15) 

with 9i — 9a + iTr /n Sl i = {0, 1, ...,n s — 1}, and 9q parametrising the azimuthal degeneracy 
of solutions. 

The second one, characteristic of the travelling wave family studied here, is a combined 
shift-reflect symmetry. Solutions invariant under this symmetry operation, 

Su(x) — S(u, v, w)(r, 9i + 9, z; t) — (u, —v, w)(r, 9i — 9, z + ir/n; t), (2-16) 

are left unaltered when shifted half a wavelength downstream and then reflected with 
respect to any of n s diametral planes tilted with 9{. 



A number of highly symmetric solutions possessing both symmetries (Pringle et al. 



2009) have been identified and continued to very low Re numbers. Some of the less 
symmetric waves are known to emerge from symmetry breaking bifurcations of these 
highly symmetric waves. 

All symmetries can be enforced upon computation. This obviously applies to time 
evolution, computation of relative equilibria using Newton iteration and stability analysis 
of relative equilibria via Arnoldi iteration. Nevertheless, we have chosen to only restrict 
the azimuthal wave number n s , apart from the invariance due to axial periodicity, to avoid 
inevitable instability to azimuthally non-symmetric flows. Every additional symmetry has 
been exposed to eventual break-up. 

2.3. Computation of relative equilibria, continuation and stability analysis 

The spectral method discussed in section |2.1| is extremely versatile and very well suited 
for straightforward coding of a Newton method for the computation of relative equilibria, 
for their continuation in parameter space and for linear stability analysis using Arnoldi 
iteration. Some of these usages of the spectral method were briefly described in |Meseguer| 



et al. (2007) but will be developed further here. 
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2.3.1. Computation of travelling/rotating waves 

The spectral representation of a generic travelling and/or rotating wave reduces to 

a(t) = a tw e- i(K ^ + " 3 ™ ce)t , (2.17) 

where a tw is the state vector corresponding to the time-independent structure of the 
wave and c z and eg are the axial and azimuthal advection speeds, respectively. The 
travelling-waves family studied here has no rotation and, therefore, eg — 0. 

Substitution of ( |2.17[ ) into ( |2.7|2.8[ ) yields a non-linear system of algebraic equations 

F(a tw ,c z ,ce, f) = 

Q ■ o,tw = o , 1 , 

^(«tw) = ^(a°) ' [Z - W) 
(j)g(a tw ) = cj)g(a°) 

where the last two scalar equations impose the phase on any two suitable coefficients so 
as to lift the translational and rotational degeneracy of the solution. The vector function 
F is given by 



F(a tw ,c z ,c e , f) = ( — B 1 + B 2 + i(nlc z + n s nc g )A ) a tw - 6(a tw ,a tw ) + f C. (2.19) 



System (2.181 can be solved via Newton iteration. At step k, the linear system to be 



solved takes the form 



D a G k 

Q T 



D Cz G 






D ce G K 






ReB^ 






c 



( Aa k 

Ac k z 



( Gk \ 




V o ) 



(2.20) 



where Stokes preconditioning (multiplication by ReB 1 from the left) has been ap- 
plied to the first row to assist convergence (Mamun & Tuckerman 19951). Accordingly, 
_ <„ _ _ 1 i i i _i . ■ — ii — - — - — ■ _ 



G h 



ReB 1 1 F(a k ,c z , Cg, f k ). The terms appearing on the first row of the left hand side 



correspond to the Jacobian of G, noted D X G, evaluated at step k: 

D a G k = l + ReB^ 1 B 2 +iRe{nlc k z +n s nc k e )B^ 1 A-ReBi 1 D a b(a k ) 



D C G K 
D c ~ e G k 



inlReB^A^ 
in s nReB^ 1 Aa k 



(2.21) 



where D a b(a k )Aa k , the linearised nonlinear term, is computed in physical space as 
(Au k • V)ti tw fc + (u tw k • V)Au fc and then projected back to spectral space (Meseguer 



& Mellibovsky 2007). Arrays <p z and (f) s correspond to the linearisation of the phase 



prescription of any two coefficients, as long as they are associated with basis elements 
that entrain axial and azimuthal structure, respectively. In the notation of Meseguer & 



Mellibovsky| (|2007[) , we have chosen to fix the phase of a^o an d a oio to u ^ the axial and 



,(2) 



azimuthal degeneracies of the solutions, respectively. The incremental vector the Newton 
method requires solving for at each step is defined as Ax k — x k+1 — x and provides the 
new guess for step k + 1. 

The solution of the linear system at every Newton step needs only be approximated by 



using inexact Krylov iteration. Generalised minimal residuals (see Fraysse, Giraud, Grat 



ton & Langou 20031 for the implementation of GMRES we use) and stabilised biconjugate 



gradient methods (BiCGSTab) can be easily implemented (see e.g. Quarternoni, Sacco 
fc Saleri|2007 and references therein) . Convergence to the exact solution of the nonlinear 



equations to numerical accuracy is guaranteed provided that a sufficiently good initial 
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guess (a , c°, c§, /°) is fed into the iteration. Convergence can be enhanced by appropri- 
ate damping (Quarternoni et al. 2007[ ) or trust-region techniques (Dennis & Schnabel 
1996). 



2.3.2. Continuation of solutions 

It is fairly easy to extend the Newton method just described to the continuation in 



Re of any branch of solutions. The explicit dependence of system (2.181 on Re allows 



this parameter to be treated as an additional unknown. Orthogonality to the Jacobian at 
some prescribed distance (pseudo-arclength) from the departing solution provides the ad- 
ditional equation needed ( Sanchez, Marques k, Lopez|2002[ ). This continuation technique 
is capable of following branches around turning points, which makes it very convenient 
in the presence of saddle-node bifurcations. Continuation in the other parameter k is 
not as straight forward. The complex dependency of the linear contribution of F on k 
requires recomputing the linear operators at every single Krylov iteration, which renders 
the method very costly. We have chosen to continue it manually by simply using the 
Newton method with discrete jumps in k, small enough so that the converged solution 
at the previous k constitutes a sufficiently good guess. This procedure cannot go around 
turning points, but combined with pseudo-arclength continuation in Re, the full family 
of travelling waves can be computed. 

2.3.3. Stability analysis 

Special care must be taken when analysing the linear stability of rotating/travelling 
waves. Waves of this type have two equivalent representations in the constant pressure 
gradient and in the constant flow-rate frames. This is so because their solid-body mo- 
tion neither modifies the mass-flux nor the mean axial pressure gradient. However, their 
destabilisation does not necessarily preserve these two global quantities simultaneously, 
meaning that linear stability will explicitly depend on which of the two is prescribed 
as constant. Since we have chosen to study pipe flow at a constant flow rate, we must 
accordingly allow for pressure variation. Let us consider an infinitesimal perturbation 
(v a ,Vf) on a given travelling/rotating wave (atw,/tw) in a reference frame moving with 
the wave: 

(2.22) 



v a c- 



a = [a tv 

f = /tw + VfC- 

The factor e A * is intended to capture the linear evolution of the perturbation in the 
vicinity of the relative equilibrium. In this sense, ( v a ,Vf ) is an eigenvector of eigenvalue 
A, which in general can be complex. Introducing (2.22) into ( 2.7|2.8 | yields the evolu- 
tion equation for the perturbation, which can be linearised around the wave solution by 
discarding all second order terms: 



B 2 +i{Klc z + n s nc g )A- D a b(a tw ))v a + Cvf = XAv c 

Q-v a = 



(2.23) 



System (|2.23|) constitutes a constrained generalised eigenproblem which can be expressed 

Ax = XMx, 



in matrix form as 
where 



A 



■j^B 1 + B 2 + i(nlc z + n s nc g )A - D a 6(a tw ) C 



A " 



M = 



T 



Q T 

and 



v a 
v f 



(2.24) 



(2.25) 
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Since we are only interested in leading eigenvalues (those with larger real part) while 
usual iterative methods are good at capturing dominant eigenvalues (those with larger 



modulus), we solve (2.24 1 in shift-invert mode 



{A-aM) 1 Mx = isx, 



A 



(2.26) 



with a real shift a applied to avoid the neutrally stable double eigenvalue at the origin 
corresponding to infinitesimal rotations and translations. 

The method used is Arnoldi iteration with Krylov-based iteration for the solution of 
the linear system involved at each step (Freitag 2007). We use the robust routines by 



Lehoucq & Scott (1996). Stokes preconditioning is employed as for the Newton method 



to speed up convergence. 



2.4. Normal form reduction 

We will be concerned here with the stability and bifurcation of travelling-wave-type solu- 
tions. Strictly speaking, a Poincare section should be defined and the bifurcation analysis 
carried out for the map defined by the resulting Poincare application. In this setting, trav- 
elling waves would be seen as fixed points and modulated waves as periodic orbits. It is 
common practice, though, to approximate the map to any desired accuracy with a system 
of differential equations (see Guckenheimer & Holmes 1983 chapter 4). Bifurcations of 
maps can thus be studied by constructing continuous-time normal forms. This comes out 
naturally when studying the stability of relative equilibria in a comoving reference frame 
(2.22). In the close neighbourhood of the relative (or degenerate) equilibrium, such an 
analysis is perfectly accurate ( Rand|19 82; Kr upa|1990 ). These travelling waves may typ- 
ically undergo Hopf bifurcations. The bifurcation is considered a Hopf bifurcation from 
a relative equilibrium rather than a Neimark-Sacker (Hopf bifurcation for maps) because 
of the special type of time-dependence involved. The bifurcating quasi-periodic solution 
is a modulated travelling wave which is just a periodic orbit in the comoving frame in 
which the travelling wave is stationary. Modulated waves that bifurcate via Hopf bifurca- 
tions from drifting waves in systems with 5*0(2) symmetry were classified in terms of the 
spatial symmetry of the drifting wave by Rand (1982). The classification was completed 
by Golubitsky, LeBlanc & Melbourne ( |2000 ) to include the spatiotemporal symmetry 
of the modulated wave. The appropriate way to investigate these solutions is to study 
the dynamics normal to the group orbit (Krupa 1990), which in our case is equivalent 
to considering the dynamics in a moving reference frame in which the travelling wave is 
stationary. 

A problem arises when, for given values of the parameters, coexistence of several waves 
occurs and, possibly, of one or more modulated waves. In this case, each wave will gen- 
erally drift with its own phase speed, so that there is no unique moving reference frame 
in which their degeneracy could be lifted simultaneously. This is easily overcome by stat- 
ing the problem back with the use of Poincare sections, combined with continuous-time 
interpolation of the corresponding maps. This of course only solves the problem as long 
as the limits of validity of the interpolation are not infringed. In particular, special care 
should be taken when studying the dynamics in the vicinity of homoclinic orbits, where 
arbitrarily small perturbations can lead to qualitatively different behaviour due to man- 
ifolds entanglement. In the absence of translational symmetry, complex dynamics would 
arise, i.e. horseshoes ( Smale|1967 Guckenheimer fc Holmes|1983 ). However, it is possible 
to show that, due to the translational symmetry of the pipe, if a point lies on both the 
stable and the unstable manifold of a saddle such as a lower-branch travelling wave, then 
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the unstable manifold is contained within the stable manifold and complex dynamics 
does not occur ( jChossat fc LauterbachpOOO ) . 

It will be argued that, in the problem at hand, there are only two active modes in the 
range of parameters studied, all other modes being slaved to these two via the centre 
manifold theorem. This, added to the translational symmetry of the pipe, results in the 
ability of capturing the dynamics within this region of parameter space with a pair of 
amplitude equations: 

x = Fi(x,y), , 22? s 
V = F 2 (x,y), 

where x and y are real amplitudes related to velocity fields in (2.1 1. The simplest expan- 
sion for Fi and F2 that is compatible with the dynamics observed upon integration of 
(2.7) is the subject of normal form theory and will be developed along the analysis. 



3. Results 

3.1. Travelling waves 

The special stability properties of the specific family of travelling waves in the 2-fold 
periodic subspace singles them out for the kind of analysis performed in this study. These 
waves, whose existence relies on the Self- Sustaining Pr ocess (SSP) (|Wa leffe 1995 1997), 
were first computed by Faisst & Eckhardt (2003) and Wedin & Kerswell (2004) using 
volume forcing homotopy. They are characterised by discrete 2-fold azimuthal periodicity 
(Z2, i.e. n s — 2) and a combined shift-reflect symmetry (S). 

For a given k and above a certain critical Re, they appear in pairs, one that is closer 
to the laminar solution, the other being further away in all respects and whatever the 
property examined (friction factor, energy contents, distance in phase space, etc). They 
are accordingly tagged as a lower-branch and an upper-branch solution, respectively. 
Figure [I] (see online movies) shows axial velocity contours together with in-plane velocity 
vectors on four cross-sections spanning half a wavelength for k = 1.52, Re = 1600 
upper-branch (top row) and lower-branch (bottom row) travelling waves. The discrete 
2-fold azimuthal periodicity can be seen in all cross-sections. Also the combined shift- 
reflect symmetry can be grasped by comparing the z = and z — A/2 cross-sections of 
either travelling wave. The reflection plane is tilted with 9q — 7r/4. The complementary 
half-wavelength is not shown, as it can be constructed by mere reflection of the given 
cross-sections. The travelling waves shown are representative of the full family in that 
they are characterised by wobbling low-speed streaks (light shading) in the pipe central 
region and fairly streamwise-independent high-speed streaks in the wall region (dark). 
Leaving aside weak secondary vortices, a couple of strong vortex pairs can be identified 
from the velocity vectors. In each of the pairs, the vortices are not of equal intensity, but 
alternate in strength along the axial direction. Except for the higher axial velocities and 
the slightly wider low speed streaks, not much distinguishes the upper-branch solution 
from the lower-branch solution. 

Radial-velocity contour-lines on a radial plane at r = 0.65 have been plotted in figure[2] 
for both travelling waves. The difference between upper-branch and lower-branch waves 
becomes clearer, the former involving much larger radial transport. These sections display 
clearly both the azimuthal periodicity and the shift-reflect symmetry of the waves. Struc- 
tures resembling these and other travelling waves have been observed experimentally in 



developed turbulence (Hof et al. 20041 as well as in numerical simulation (|Kerswell & 
Totty|[2007l |Schneider, Eckhardt \ Vollmer|[2007| |Willis fc KersweIIl[2008| . 



They appear in saddle-node bifurcations starting from Re ~ 1358.5 for the optimal 
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Figure 1. Cross-sections spanning half a wavelength (A/2 = 7t/k) of the upper-branch (top 
row) and the lower-branch (bottom row) travelling waves for n — 1.52, Re — 1600. Axial 
perturbation velocity contour-lines are equispaced with intervals of Au z = 0.1(7. Regions where 
the axial velocity is faster or slower than the laminar parabolic profile are coded in dark and 
light, respectively, with the zero contour-line marked with a dashed line. In-plane velocity vectors 
are also shown. 
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Figure 2. Radial sections at r = 0.65 of (a) the upper-branch and (b) the lower-branch travelling 
waves for k, — 1.52, i?e = 1600. Radial velocity contour-lines are equispaced with intervals of 
Au r ~ 0.007217. Regions where the flow goes towards the walls or the pipe center are shaded 
dark and light, respectively. 



wave number k ~ 1.55 ( Wedin fc Kerswell|2004 ). Lower-branch solutions exhibit a single 
unstable eigenmode when considered in the azimuthal subspace where they live, which 
makes them accessible through time integration if appropriately combined with edge 
tracking techniques. In this sense, they are edge states (Schneider et al. 2007) of the 
2-fold azimuthally-periodic pipe, and they constitute an example of non-uniqueness of 
the edge state, since they coexist as attractors within the critical threshold, and for some 
range of n and Re, with another family of travelling waves that are both shift- reflect and 
mirror symmetric and that have also a single unstable direction ( Duguet et aZ.|2008 ). The 
stable manifold of these waves separates, albeit only locally, initial conditions that decay 
uneventfully from others that undergo turbulent transients. For this reason, lower-branch 
solutions are believed to play a central role in the transition process. 

As already pointed out, lower-branch solutions, when continued down in Re, merge 
with upper-branch solutions in a saddle-node-type bifurcation. Around the saddle-node 
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Figure 3. Rightmost eigenvalues of (a) n — 1.50 and (6) k — 1.70 travelling waves at their 
respective saddle- node bifurcation (SN). The direction of the zero-crossing from lower to upper 
branch is indicated with an arrow. The eigenvalues marked as filled circles have been followed 
across the saddle-node point and beyond for both (c) k = 1.50 and (d) k = 1.70. Saddle-node 
(SN), node-focus (NF) and Hopf (H) points are identified with circles. 



point, a real eigenvalue must cross the imaginary axis, either adding or subtracting an 
unstable direction. As a matter of fact, both scenarios occur and the stability properties 
of upper-branch solutions depend on the axial wavenumber k. Upper-branch long waves 
are unstable, whilst short waves are stable. The spectrum of a long (k — 1.50) and of a 
short (n = 1.70) wave are shown in figures^a) and[3]j6), respectively, at their respective 
saddle-node points. The arrow below the bifurcating eigenvalue points in the crossing 
direction when followed from lower to upper branch. For k = 1.50, the solution goes from 
one to two unstable directions, while for n = 1.70 it goes from one to nil. Figures [3]( c) 
and[3j<i) show the trajectories followed by the two rightmost eigenvalues (filled circles 
in figures [3)3,6) when continued across the saddle-node (SN) and beyond, for n = 1.50 
and k = 1.70, respectively. In both cases, the eigenvalues collide and go complex in a 
node-focus (NF) transition, but they can do so at either side of the imaginary axis. When 
followed further to higher Re along the upper branch, long and short solutions undergo 
a stabilising or destabilising Hopf bifurcation (H) , respectively, implying the existence of 
modulated waves (relative periodic orbits). 

Cross-sections of both travelling waves at their respective saddle- node points have been 
plotted in figures 4|a) andQc), together with their corresponding crossing eigenfunctions 
(figures |4^) and[4i). It is clear from the contour plots that the effect of varying n on the 
solutions does not affect the flow structure dramatically, aside from the stretching and 
the stability properties. It is also remarkable how the eigenfunctions corresponding to 
two different eigenvalues (figures |4]c,rf) look so similar to each other, even for different n. 
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Figure 4. (a) k = 1.50 and (c) n — 1.70 travelling waves at their respective saddle-node bifur- 
cation points (SN in figure Eb. Also shown are the bifurcating eigenfunctions for (b) k = 1.50 
and (d) k = 1.70. Axial perturbation velocity contour-lines are equispaced in intervals of 
Au z = 0.1(7. 




This is a signature of the tendency to become parallel, which ultimately leads to collision 
in a node-focus (NF) transition. Both eigenfunctions preserve the same symmetries as 
the travelling waves, so that all solutions arising from their bifurcation must necessarily 
possess the same symmetries. 

Wavenumber continuity suggests that there might be an intermediate value of k for 
which the node-focus point can be brought to the origin of the complex plane. Appro- 
priate tuning of k and Re would therefore help unfold a codimension-two double zero 
bifurcation, a complete study of which is due to Takens (1974) and to |Bogdanov (1975). 
As a matter of fact, the trajectories of the complex eigenvalues, following a left opening 
parabola plus an eventual right turn, indicate that for a certain range of k an additional 
pair of crossings will take place, thus implying additional Hopf bifurcations closely related 
to the appearance of periodic orbits in a double-zero bifurcation. 



3.2. Modulated travelling waves 

Figure [5j a) depicts the trajectories of the relevant eigenvalues for k = 1.52. They are 
extremely close to those observed for k = 1.50, except that the nose is getting tangent to 
the imaginary axis. The implications of this will be discussed later on. Let us focus on the 
Hopf bifurcation (H) for the time being. The crossing takes place at Re = 1775.4 on the 
upper branch, with the arrows pointing towards higher Re. The real and imaginary parts 
of the bifurcating eigenfunction are shown in figure [5^6) only for completeness, since it 
is difficult to extract from them relevant information on the nature of the instability 



Takens-Bogdanov of travelling waves 13 
(a) (b) 




Figure 6. (a) Fourier transform of the non-axisymmetric streamwise-dependent modal energy 
contents (e3d) time-stamp for the (n,Re) = (1.52,1600) modulated wave. The time signal 
is plotted in the inset, (b) Advection speed (c z , solid line) and mean axial pressure gradient 
normalised by its laminar value ((Vp) z , dashed line) along a full period. The values for the 
(ft, Re) — (1.52, 1600) upper-branch travelling wave are indicated by the horizontal line. 



that they carry. As already noted, the instability does not break any of the symmetries 
exhibited by the bifurcating travelling wave. 

The existence of a Hopf bifurcation implies the appearance of a time-periodic solution, 
in this case a modulated wave or relative periodic orbit. Since the crossing is leftwards as 
Re is increased, it necessarily carries with it a stabilisation of the upper-branch travelling 
wave, which was unstable at onset. The Hopf bifurcation happens to be supercritical in 
the sense that the first Lyapunov coefficient is negative ( Kuznetsov||1995 ). This we infer 
from the fact that time evolution unveils a stable branch of modulated waves pointing 
towards decreasing Re. 

Simple time evolution starting from the (k, Re) = (1.52, 1600) unstable upper-branch 
travelling wave (figure [I] top row) departs in an oscillatory fashion, as predicted by the 
complex pair of unstable eigenvalues, to nonlinearly converge onto a modulated wave. 

Travelling waves, which are degenerate or relative equilibria, and solutions bifurcating 
from them have a pure frequency associated to the advection speed (solid-body trans- 
lation). Local quantities such as point velocities or pressures necessarily reflect on this 
frequency. However, since it reflects a neutral direction, it can be suppressed by restating 
the problem in a comoving frame. Global quantities such as modal energies or volume 
averaged fields naturally overlook solid-body rotation and translation, making them suit- 
able for a decoupled analysis in the direction orthogonal to the degenerate drift. Thus, 
as justified before, the bifurcation analysis of travelling waves can be carried out analo- 
gously to that of fixed points, as long as special care is taken in the neighbourhood of 
homoclinic connections ( Rand|1982 Golubitsky et aL|2000 |. The Hopf bifurcation adds a 
modulational frequency to the pure translational frequency. As a result, global quantities 
cease to be constant and oscillate with this frequency. The Fourier transform of the en- 
ergy contained in non-axisymmetric streamwise-dependent modes (£3d) has been plotted 
in figure^ a), with the time signal shown in the inset frame. The spectrum reveals that 
the solution has a strong mean component and a peak angular frequency at loq ~ 0.0219, 
corresponding to a period T = 2tt/u> = 286.7. This period is extremely long when com- 
pared with the streamwise advection time-scale, which is of order 2tt/{kc z ) ~ 3. The 
signal is not strictly sinusoidal and some energy is spread among a number of harmonics 
of wo- A full period of the axial advection speed (c z ) and of the mean axial pressure 
gradient normalised by its laminar value ((Vp) 2 ) has been plotted in figure If 6). Both 
are global quantities and, as expected, oscillate with the modulational frequency cj . The 
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Figure 7. From left to right, t = 50, 125, 175 and 225 snapshots of the ((«, 7?e) = (1.52, 1600) 
modulated wave, as indicated by the vertical dotted lines in figure [6j b). Shown are the z = 
(top row) and z = A/4 (bottom row) cross-sections with axial velocity contour-lines equispaced 
in intervals of Au z = 0.1U. To avoid drift due to streamwise advection, snapshots are taken in 
a comoving frame travelling at the instantaneous advection speed from figure j6[ b). 



horizontal line indicates the constant value of these quantities for the unstable upper- 
branch travelling wave from which the modulated wave bifurcates. The modulated wave 
orbits in phase space around the travelling wave, so that the instantaneous values of 
the former oscillate around those of the latter. The differences between time-averaged 
quantities and travelling wave values are purely due to nonlinear effects, since for this 
Re and k we are already some distance away from the linear influence of the instability. 

Figure [7] (see online movie) depicts a couple of cross-sections [z = 0, top row; z = A/i, 
bottom row) of the modulated wave at four different time instants along one period (from 
left to right, t = 50, 125, 175 and 225). The modulation is obvious from the series of 
snapshots yet not very prominent, as could be expected from the mild oscillation (within 
±5%) of the time series in figure's). Comparison with figure[l](top row) reveals that the 
overall flow structure of the travelling wave is preserved on average but is now dependent 
on time. 

Radial sections showing radial velocity contours have been plotted in figure [8] (see 
online movie) to convey in a clear way the modulation of the wave along a full period. 
It also confirms that all symmetries are preserved. It will soon become clear that the 
snapshot at t — 50 is about the closest approach to the lower-branch travelling wave, 
while all of them are approximately equally distant from the upper-branch wave. 

Straightforward i?e-continuation of the n — 1.52 modulated waves by time evolution 
reveals that they bifurcate off the travelling waves upper branch precisely at the Hopf 
bifurcation point at Re = 1775.4 (H in figure [5)3). Figure [9] shows phase map represen- 
tations of the modulated waves for a discrete set of Re along the branch. Two global 
quantities, £30 and (Vp) z , have been chosen to get rid of axial drift effects. Increasing 
Re from Re = 1600 (towards lighter gray), the modulated wave (solid closed loops) 
shrinks and closes on the upper-branch travelling wave (filled circles), to end up merging 
at the Hopf point (Re — 1775.4, not shown). In the opposite direction, the wave also 
starts shrinking but the process stops, slightly reverses and track of the wave is lost, as 



Tokens- Bogdanov of travelling waves 
t = 50 t = 125 

H 



15 




Figure 8. Radial sections at r — 0.65 of the (k, Re) = (1.52, 1600) modulated wave. Snapshots 
are taken in a reference frame moving at the wave speed. Radial velocity contour-lines are 
equispaced with intervals of Au r = 0.0072(7. Regions where the flow goes towards the walls or 
the pipe center are shaded dark and light, respectively. 
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Figure 9. Phase map projections of travelling and modulated waves along a .Re-continuation 
for k = 1.52. The axes represent non-axisymmetric streamwise-dependent energy (esd) against 
mean axial pressure gradient ((Vp) z ). Different gray intensities correspond to different Re as 
explained in the legend. Lower- and upper-branch travelling waves are marked with pluses and 
filled circles, respectively; modulated waves with solid lines and escaping trajectories with dotted 
lines. 



pointed out by the diverging dotted trajectory. This trajectory wanders around a ghost 
of the disappeared modulated wave for some time, then spirals away until it comes close 
to the lower-branch travelling wave (plus signs) and leaves the region towards the lam- 
inar flow. Close analysis seems to indicate that the branch bends and turns around a 
fold (saddle-node) of cycles picking up an unstable direction. The phase map picture 
is, in this sense, incomplete. For each of the stable cycles (nodal branch) there must 
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Figure 10. (a) Eigenvalues trajectories for ft = 1.53. (6) Phase map projections of travelling 
and modulated waves along a .Re-continuation for n — 1.53. See the caption of figure [9] for 
explanations. 



be an unstable cycle (saddle branch). As Re increases from the fold-of-cycles point, the 
nodal branch approaches the upper branch of travelling waves until eventual collision at 
the Hopf point, while the saddle branch grows larger and larger. It will be argued later 
that, unless some exotic symmetry-breaking event takes place, the unstable modulated 
wave must end up colliding with the lower branch of travelling waves in a saddle-loop 
bifurcation, where it ceases to exist without any other consequence. 

As already mentioned, decreasing Re causes both branches to merge in a fold-of-cycles, 
leaving unstable travelling waves as the only existing solutions. Any trajectory issued 
from the vicinity of the upper-branch wave spirals away, temporarily orbits around the 
ghost of the disappeared modulated wave, and then continues spiraling away until the 
trajectory manages to go around the lower-branch travelling wave and finally leaves 
towards the basic laminar profile (ssdj (Vp)*) = (0, 1) along the only unstable direction 
of the travelling wave. 

Continuation in the geometrical parameter k causes the node-focus (NF) transition 
point to move around without noticeable impact on the eigenvalues i?e-trajectories, as can 
be concluded from figures [3]( c, (i) and[5]ja). It becomes clear that for wavenumber slightly 
beyond n > 1.52 the trajectories will experience an additional couple of zero-crossings, 
adding a pair of Hopf bifurcation points to the already existing one. This is illustrated 
for k = 1.53 in figure [lu| a), where the successive crossings along the upper branch have 
been labelled Hi, H2 and H3 as Re is increased. Continuation of the modulated waves 
branch for k = 1.53, shown in figure [To^ 6) , reveals that they connect supercritically to 
the upper branch of travelling waves in both directions at the two Hopf points that are 



at a larger distance from the saddle-node bifurcation, i.e. H 2 and H 3 of figure 10 [a). The 
fold-of-cycles that was present at k = 1.52 is no longer present at k = 1.53. 

We shall extend this notation, Hi, H2 and H3, to denote the three potential Hopf 
crossings for any k in the range studied. It is clear that for decreasing n, Hi and H 2 
approach, collide and disappear leaving H3 alone. Instead, for increasing k it is H 2 and 
H3 that approach and collide, leaving Hi alone as long as the node-focus point remains 
on the right side of the complex plane. 

Adding the i?e-dependence helps to clarify the situation. This is done in figure [XT] 
where £ 3 d has been plotted as a function of Re and k (varying gray intensity as indicated 
by the legend) for both travelling waves (solid, dotted and dashed thick lines for 0, 1 
and 2-dimensional unstable manifolds) and modulated waves (filled circles indicating 
minimum and maximum values of £30 along a period). Lower-branch travelling waves 
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Figure 11. Non-axisymmetric streamwise-dependent energy (e3d) as a function of Re and k 
(different gray intensities as indicated in the legend). Lines represent travelling waves: solid, 
dotted and dashed for 0, 1 and 2 unstable directions. Filled circles on the curves represent Hopf 
points. The rest of circles represent stable modulated waves. 



have a single unstable direction. For all cases shown, as they turn around the saddle-node 
point, a second unstable direction appears. Low-k (e.g. n — 1.52) waves remain unstable 
until they regain stability at Hopf bifurcations, whence supercritical branches of stable 
modulated waves are issued, as already discussed for figure [9} Traces of the turning point 
that we identified with a fold-of-cycles are clear from the bending tips of the modulated 
wave branches. High-K (e.g. k = 1.60) waves regain stability almost immediately at 
Hopf bifurcations but no stable branch of waves could be detected. Intermediate-K (e.g. 
k = 1.53) waves exhibit a mixed behaviour by regaining stability early on the upper 
branch, but losing it again temporarily for some range of Re. The whole intervening 
unstable range is nevertheless covered by a stable branch of modulated waves connecting 
both Hopf points, as illustrated in figure |Tu| b). At even higher k (not shown), travelling 
waves become stable at the saddle-node, as was the case for k = 1.70 of figures [3^&,<i). 
Returning to upper-branch travelling waves that were unstable at onset, the branch of 
modulated waves that must be born at their early stabilisation could not be detected 
with time evolution. This is a clear sign that they may bifurcate subcritically, with the 
branch pointing in the increasing Re direction. 

Figure 12 shows the angular frequency w of the modulated waves as a function of 
Re and k. Only a subset of k has been plotted for clarity. And, of course, only as long 
as the branch is stable, i.e. accessible through time evolution. Along with the waves, 
the locus of Hopf bifurcation points (thin dotted line with crosses) has been added by 
plotting the imaginary part of the bifurcating eigenvalues. Along this curve, n varies 
continuously but not monotonically. The curve has been split in three monotonical seg- 
ments corresponding to the up to three independent zero-crossings that the complex pair 
of eigenvalues might undergo (Hi, H2 and H3). Naturally, the modulational frequency 
of the bifurcating modulated waves is dictated by the imaginary part of the eigenval- 
ues at the Hopf bifurcation. This figure strongly endorses the fold-of-cycles hypothesis 
for the non-reconnecting branches. The decreasing frequency around the fold-of-cycles is 
also suggestive of the eventual saddle-loop collision with lower-branch travelling waves, 
which demands a vanishing frequency (period going to infinity) as the modulated wave 
transforms into a non-robust homoclinic loop and disappears. 
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Figure 12. Angular frequency (uu) of the modulated waves and the Hopf points (filled circles). 
The thin line with crosses represents the imaginary part of the bifurcating eigenvalues along the 
locus of Hopf bifurcation points. The three segments, labelled Hi, H2 and H3, correspond to 
intersections of the Hopf curve with iso-K lines. 



It is clear that for 1.52 < n < 1.53 a change of behaviour takes place. Modulated wave 
branches go from reconnecting to the Hopf curve to bending and losing stability. It seems 
therefore reasonable that for some intermediate value, a branch of solutions may become 
tangent to the Hopf curve, precisely at the k value for which the couple of additional 
eigenvalue crossings start occurring, i.e. where Hi and H2 merge. Slightly above this 
value, the curve would split in two. One of the segments connects supercritically two 
of the three existing Hopf points (H2 and H3). The other segment would be issued 
from the third Hopf point (Hi) and, for purely geometrical reasons, it would have to 
do so supercritically towards lower Re. Unfortunately, the room left in parameter space 
for the stable branch of waves before it bends back and becomes unstable is tiny and 
very close to the Hopf bifurcation line where the relevant complex pair of eigenvalues 
exhibits a vanishing real part. As a consequence, transients are far too long to allow 
computation of the stable modulated waves in a feasible time-scale. This ultimately 
renders the numerical determination of the borders of existence of such solutions within 
this region extremely difficult. For even higher k, but not exceeding the value for which 
upper-branch travelling waves are stable at onset, modulated waves can only bifurcate 
subcritically from Hi towards increasing Re and with vanishing frequencies, as required 
by the double-zero bifurcation scenario that takes place. 

3.3. Unfolding of the Takens-Bogdanov bifurcation 

The 2-fold azimuthally-periodic, shift-reflect symmetric family of travelling waves has 
been continued both in Re and K, and systematic linear stability analysis has been carried 
out. In this manner, the locus of the saddle-node bifurcation between upper- and lower- 
branch travelling waves, and that of the Hopf bifurcation of upper-branch travelling waves 
has been tracked in parameter space. Results are summarised in figure [13] To the left of 
the saddle- node bifurcation lines (SN+ and SN_), there exist no travelling waves. These 
appear in saddle-node bifurcations to the right of the saddle-node lines. Upper-branch 
waves produced in SN + are stable nodes, whereas those appearing in SN_ are unstable, 
all lower-branch travelling waves being unstable with a single unstable direction. Almost 
immediately after springing into existence, upper-branch travelling waves experience a 
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Figure 13. Bifurcation curves in (k, _Re)-parameter space. SN+ (solid line with crosses) and 
SN_ (dashed line with crosses) represent stable and unstable saddle-node bifurcations of trav- 
elling waves. H (solid line with circles) is the Hopf bifurcation curve of upper-branch travelling 
waves. TB is the Takens-Bogdanov (double-zero) bifurcation point. Stable modulated waves are 
indicated with triangles, with FC (dotted line) representing a fold-of-cycles. The inset shows a 
transformed zoom of the dashed region, with new coordinate axes ktb — k and ResN(n) — Re- 

node-focus transition (NF of figure |3| . This non- bifurcation curve is not shown in the 
figure, as it is extremely close to the saddle-node curve. SN + , SN_ and NF emanate 
from a Takens-Bogdanov or double-zero bifurcation point (TB). At this point, two real 
eigenvalues collide exactly at the origin of the complex plane and become a conjugate 
pair. From this very same point, a Hopf bifurcation line (H) is issued downwards to 
the right of SN_. Initially, the lower is k, the further away H gets from SN_. This is 
a consequence of the left-opening parabolic shape of the complex eigenvalues (figure [3]) . 
As k is reduced, NF moves away from the origin and the distance along the trajectory, 
which is related to Re, before crossing gets longer. Further along H, a couple of successive 
upward and downward bendings occur. The three segments into which they split the curve 
can be identified with Hi, H2 and H3, the three potential crossings the eigenvalues may 
undergo at any given k. Above this line, upper-branch travelling waves are stable. Below 
it they are always unstable. As a consequence of its last downwards turn, for sufficiently 
large Re all upper-branch solutions are left stable as far as this bifurcation diagram is 
concerned and no secondary instabilities take place. To complete the phenomenological 
description, modulated waves have also been added into the picture (empty triangles). 
They occupy a region enclosed between the Hopf line and another line defined as the locus 
of turning points where modulated waves cease to exist. This line is a bifurcation curve 
corresponding to a fold-of-cycles (FC), whereby stable modulated waves turn around and 
become unstable. 

4. Discussion 

It is at this point convenient to define an alternative couple of parameters: «tb — k and 
i?esN(ft) — Re, where ktb is the value of k at the Takens-Bogdanov point and Reg^(K) is 
the value of Re at the saddle-node point for each given k. The transformation has been 
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applied to the region indicated by the dashed square in figure p~3| and plotted in the inset. 
The resulting bifurcation diagram is reminiscent of the double-zero bifurcation described 
in Wiggins (2003). Of the two equivalent normal forms proposed for this bifurcation by 



Takens ( 1974 ) and Bogdanov ( 1975 ), the former can be applied to the present bifurcation 
by simply taking \x\ ~ ktb — k and ji2 ~ i?esN(«) — Re in Wiggins (2003[): 



x = y, 



M22/- 



bxy, 



(4.1) 



with b = +1. A detailed analysis of the alternative normal form can be found in |Kuznetsov| 



(1995). 



This second-order normal form provides an accurate description of all the phenomenol- 
ogy observed in the immediate vicinity of the Takcns-Bogdanov point and gives addi- 
tional information not accessible with the methodology here deployed. From the Takens- 



Bogdanov (TB) point (inset of figure 13), stable (SN+) and unstable (SN_) saddle-node 



lines extend along the vertical axis downwards and upwards, respectively. Issued from 
the same point is a subcritical Hopf curve (H) with a left-opening half-parabolic shape, 
to the left of the unstable saddle-node line. The unstable node created at the saddle- 
node immediately becomes a focus and later recovers stability upon crossing the Hopf 
curve leftwards. An unstable branch of unstable periodic orbits originates surrounding 
the stabilised equilibrium. Close analysis of the bifurcation reveals that the unstable pe- 
riodic orbits grow large and collide with the unstable saddle in a saddle-loop homoclinic 
bifurcation, not shown in the figure. The unstable periodic orbits disappear with no con- 
sequence for the saddle equilibrium, leaving it, together with the stable focus, as the only 
surviving solutions. A saddle-loop curve, also branching off the Takens-Bogdanov point, 
must therefore be present, although it cannot be detected with the methods available. 

Analysis of the second-order normal form is nevertheless insufficient to account for the 
existence of stable modulated waves and the back and forth bending of the Hopf curve. 
It is however possible to analyse in some detail what the effect of higher-order terms 
would be if the conditions for a non-degenerate double-zero bifurcation were nearly vio- 
lated and higher-order terms were not negligible some distance away from the bifurcation 
point. A complete analysis of the codimension-3 bifurcation resulting from a vanishing 



second-order coefficient (6 in 4.1 ) can be found in Dumortier et al. ( 1991 ). It is unlikely, 
unless some symmetry at play cancels exactly some second-order coefficients, which is 
not the case (the eigenfunctions involved preserve all symmetries of the solution), that 
simple tuning of two parameters will cross a codimension-3 point. Therefore, it is not 
the codimension-3 point we are interested in, but a perturbation around it that ren- 
ders higher-order terms important at some distance from the codimension-2 point under 
analysis. 



Such an analysis was carried out by Barkley ( 1990 ) when studying similar phenomena 



in plane Poiscuille flow, i.e. the appearance of modulated waves in channel flow. Previous 
studies ( Pugh fc Saffman||1988 Soibelman &i Meiron||1991 1 had reported the onset of a 
modulational instability of travelling waves when the flow was driven by a constant pres- 
sure gradient. No instability was detected for the same waves when the flow was driven 
at constant mass-flux. Introducing a new parameter allowing to connect continuously 
both possible boundary conditions, a bifurcating scenario arose as a natural explanation. 
It was suggested that while a double-zero bifurcation could explain the appearance of 
unstable modulated waves, the simple inclusion of third-order terms to the second-order 
normal form could explain the existence of stable modulated waves for some parameter 
range. According to normal form theory, the extended normal form needs only include 



a couple of third-order terms ( Guckenheimer & Holmes 1983 Knobloch 1986 loos & 



Adelmeyer] 1998[ ) 
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x — U, 2) 
y = pi + p 2 y + x 2 + kixy + k 2 x 3 + k 3 x 2 y, 

where appropriate scaling can be performed such that £3 = ±1. The second order co- 
efficient ki is required to be sufficiently small so that third-order terms come into play, 
but non-zero so that non-degeneracy of the regular Takens-Bogdanov is preserved for 
sufficiently small p\ and p 2 . 

We shall not reproduce here the analysis of the third-order-extended normal form for 



the Takens-Bogdanov bifurcation. Instead, we refer the reader to Barkley ( 1990 ) for a 
detailed unfolding, and pick up on the main results that are relevant to the present 
case. It can be shown that appropriate tuning of the coefficients enables this normal 
form to capture some of the phenomenology observed but so far unexplained. Namely, 
setting k 3 = +1 and ki > 0, gives the right eigenvalues path up to the first right turn 
after the node-focus transition (figures [3]c, d and[5]a). As a consequence, the Hopf curve 
experiences a first turning point as seen in figure |13| allowing for multiplicity of Hopf 
points at any given k above some minimum value. A most relevant result concerns the 
nature of the Hopf curve itself. It has already been mentioned that in the vicinity of 
the Takens-Bogdanov point, the Hopf curve happens to be subcritical in the sense that 
the first Lyapunov coefficient is positive and bifurcating periodic orbits are unstable. 



Restricting (4.2) to the Hopf curve and with an appropriate change of variables and 

(4.3) 



rescaling, we obtain 

j = p((3 1+ (3 2 p 2 + S p 4 ) + 0(p T ), 



6 = 1, 

with s = +1, which is the normal form for the degenerate Hopf or Bautin bifurcation 



(Kuznetsov 1995). (3 2 ~ k% — 3k 2 p 2 is the first Lyapunov coefficient and governs the 
change in criticality along the Hopf curve. Taking k 2 > will have f3 2 vanish for some 
value of p 2 , leaving the Hopf curve supercritical beyond some point sufficiently far from 
the Takens-Bogdanov point. This is exactly what is needed to account for supercritical 
bifurcation of stable modulated waves. Careful analysis of the Bautin bifurcation also 
requires that the stable modulated waves turn around in a fold-of-cycles (FC), whenever 
they cannot branch to another Hopf point, becoming unstable. Since the original unstable 
waves collided with the saddle in a saddle-loop bifurcation in the frame of the Takens- 
Bogdanov bifurcation, it is to be expected that they will continue exhibiting the same 
behaviour even if they now initially bifurcate supercritically. Thus, even though we could 



only follow stable modulated waves to their fold bifurcation (figure 11), we can now 
reliably comment on their fate: they should most certainly disappear in saddle-loop 
bifurcations. This was already justified upon explaining figure [T2| where the trend of the 
angular frequency suggested that the period could go to infinity after turning around, as 
the saddle-loop bifurcation requires for the homoclinic connection in which they convert 
before disappearing. The Bautin point must be somewhere along the Hopf curve before 
the first turning point. As it happens, supercriticality is only extremely mild until about 
the turning point (the locus of the fold-of-cycles appears to be very close to the Hopf 
curve), resulting in our inability to compute the stable modulated waves that must be 
present. This is so because close to the Hopf curve, the dynamics are extremely slow and 
the transients accordingly long. These are the waves that are expected to turn around, 
become unstable and are destroyed in saddle-loops as predicted by third-order terms. 
Beyond the turning point, supercriticality means branching towards lower pi (higher 
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Figure 14. Attemptive bifurcation diagram in (ktb — a, ResN(^) — i?e)-parameter space. Rep- 
resented are the two codimension-2 points TB (double zero or Takens-Bogdanov) and B (Bautin 
or degenerate Hopf). B splits the Hopf curve (H) in two segments, one subcritical (h_, dashed 
thick line), the other supercritical (h+, solid thick line). From TB, stable and unstable saddle- n- 
ode curves (SN+ and SN_), the Hopf curve and a saddle-loop curve (SL, dotted line) are born, 
together with stable and unstable node-focus transition curves (NF+ and NF_, dash-dotted 
lines). A fold of cycles (FC) is issued from B. Numbers in parenthesis indicate different regions. 



Re). These waves we have computed, but their fate crucially depends upon the second 
turning of the Hopf curve, and third-order terms tell us nothing about it. 

Adding fourth-order terms is a straightforward but tedious calculation that we shall not 
reproduce here. It cannot be expected that fourth-order amplitude equations capture the 
behaviour reliably so far from the Takens-Bogdanov point, yet there is strong numerical 
evidence that fourth-order terms have a say. They allow for an extra turning point to the 
Hopf curve, without altering any of the previous predictions. As a result, stable modulated 
waves born between the two turning points can disappear in a second supercritical Hopf 
after the second turn, as illustrated in figures [TT] and [13] Further down the Hopf curve, 
beyond the second turning point, stable modulated waves either behave us just explained 
or progress stably until they meet the fold-of-cycles (FC) and turn around. This curve 
must emerge from the Hopf curve precisely at the Bautin point, the existence of which 
can be inferred but its location not determined precisely with the tools at hand. 

Figure 14 completes the picture just sketched in an endeavouring diagram for all bi- 
furcations in parameter space. Modified parameters (ktb — Res^(n) — Re) have been 
chosen for simplicity. The full diagram is organised around the Takens-Bogdanov (TB) 
codimension-2 point. This point separates two branches of saddle-node lines (SN + and 
SN_). At SN + a stable node and a saddle equilibrium coalesce, while crossing SN_ gen- 
erates an unstable node and a saddle. The node, stable or unstable, experiences a quick 
transition and becomes a focus, stable or unstable, in a node-focus transition (NF+ and 
NF_). These are non-bifurcation curves. From TB, a Hopf bifurcation curve H springs 
quadratically to the same side as SN_. Sufficiently close to TB periodic orbits bifur- 
cate subcritically from the node at H (h_). Further away, H changes criticality at a 
codimension-2 Bautin point (B, degenerate Hopf) and becomes supercritical (h + ). At B, 
a fold-of-cycles curve (FC) is born to the h + side of H and parting from it quadratically. 
Between h + and FC, stable periodic orbits exist. At h_ and FC, unstable periodic orbits 
are generated. They must all cease to exist, at least in the vicinity of TB, in saddle-loop 
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Figure 15. Phase map diagrams corresponding to all bifurcation curves and numbered regions 
in figure [l4"| Small circles correspond to relative equilibria (filled for stable equilibria), and closed 
thick lines to relative periodic orbits (solid for stable, dashed for unstable orbits). Thin lines 
with arrows indicate stable and unstable manifolds, homoclinic and heteroclinic connections. 



bifurcations (SL). At SL, the unstable periodic orbit has grown large and collides with 
the saddle producing a non-robust homoclinic cycle. The periodic orbit vanishes while 
the saddle experiences no change whatsoever. Among all bifurcation curves drawn in 
figure |14[ SL is the most uncertain. Since it runs on the saddle, independent from the 
nodal equilibrium, there is no restriction on whether this line can cross H. Thorough 
analysis of the Takens-Bogdanov bifurcation shows that it cannot cross h_ and that it 
actually departs quadratically from SN_ , but higher-order terms can easily make it cross 
h + . Unfortunately no linear stability analysis of the saddle can cast light on the issue, 
nor can the unstable periodic orbits be tracked with time evolution. However, the picture 
just sketched would not change much. 

Figure [T5| shows phase map diagrams of the flow projected on the 2-dimensional centre 
manifold. The labels indicate to which region in figure [14] each phase map corresponds. 
Making a clockwise round-trip can help clarify the picture. Starting from (1), where no 
solutions exist, and crossing SN + into (2) a saddle and a stable node are created. The 
stable node turns into a stable focus through NF + in (3). In (4), an unstable periodic 
orbit exists, that has been created in SL. Upon crossing h + into (5) a stable cycle appears 
around the focus, which becomes unstable. The coexistence of an unstable and a stable 
cycle in (5) lasts until they collide and disappear in FC, leaving only the saddle and 
an unstable focus in (6). It is possible to arrive to (6) directly from (4). In this case, 
the unstable cycle just merges with the stable focus in h_ to leave an unstable focus. 
Crossing NF_ the unstable focus transforms into an unstable node in (7). The unstable 
focus coalesces with the saddle through SN_ back into region (1). The possible crossing 
of SL and h + would only affect the order in which the stable cycle merges with the focus 
and the unstable cycle collides with the saddle, when going from (5) to (3). 

It is at this point necessary to justify why we have not made use of the discrete 
azimuthal periodicity or of the shift-reflect symmetry in deriving the normal form. The 
fact is that we are entitled to ignore these symmetries because they play no significant 
role in the problem at hand. If any of the solutions had broken any of the symmetries, 
these would have needed to be built into the normal form. Since there is no numerical 
evidence from either time evolution or linear stability that this happens, these symmetries 
can be safely ignored. 
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5. Conclusion 

A complex bifurcating scenario in pipe Poiscuille flow, involving the appearance of 
travelling as well as modulated travelling waves, has been thoroughly analysed within 
a given azimuthal subspace, taking advantage of the special stability properties of the 
solutions it is home to. These stability properties make this subspace particularly fit for 
the search and analysis of all kinds of exotic behaviour including, but not limited to, the 
double-zero bifurcation here unfolded. This can help place a lower bound to the extreme 
complexity that must certainly be present in full space pipe flow. 

The particular family of waves studied involves the appearance in Hopf bifurcations of 
modulated waves that are either unstable at onset or become unstable in a fold-of-cycles. 
It cannot be guaranteed that the unstable modulated waves created in the fold-of-cycles 
eventually disappear in a saddle-loop bifurcation. Their ultimate fate depends on the 
absence of symmetry-breaking bifurcations that free them from the need of colliding 
with the saddle or that allow for an entanglement of the stable and unstable manifolds 
of the saddle in the vicinity of the global bifurcation. While such a behaviour has to be 
expected generically, the absence of chaotic dynamics around the saddle suggests that 
this is not the case, and that the scenario described may be fairly accurate. 

The Takens-Bogdanov bifurcation experienced by this particular family of travelling 
waves must not be considered as an exotic scenario. There is strong evidence that other 
families of travelling waves might be bifurcating in similar environments. This is the 
case of the highly-symmetric 3-fold azimuthally-periodic family of travelling waves de- 



scribed in Pringle et al. ( 2009 ) , for which the eigenvalue associated with the saddle-node 
bifurcation is known to collide with another eigenvalue and go complex immediately 
after crossing from lower to upper branch. This happens in the context of instability 
dictated by several additional unstable eigenvalues, which renders the techniques used 
in the present study unable to unfold the bifurcation. Nevertheless, a Takens-Bogdanov 
bifurcation appears to be there all the same. 

The scenario is nonetheless not general. Other families of travelling waves, such as 
the mirror-symmetric 2-fold azimuthally-periodic, do not follow this behaviour. Other 
travelling waves come from Pitchfork bifurcations from other more symmetric waves, 



e.g. the shift-reflect azimuthally non-periodic travelling wave ( |Pringle fc KerswellpOT) , 



which comes from symmetry breaking of a mirror-symmetric wave. The Takens-Bogdanov 
bifurcation simply does not apply in these cases. 

As a result of a Takens-Bogdanov bifurcation, no trace of modulated waves is expected 
to survive long enough to take part in the formation of a chaotic attractor or saddle as 
a precursor to turbulent dynamics. Nevertheless, it is interesting to point out that, far 
enough from the Takens-Bogdanov point, only stable upper-branch solutions are left. 
Thus, all the techniques here deployed are left valid for an eventual analysis of the 
bifurcation cascade they may undergo at higher Re. A detailed study of the behaviour 
of these waves at higher flow rates is currently underway. 

This work has been supported by the Spanish Ministry of Science and Technology, 
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by the Deutsche Forschungsgemeinschaft. 
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